Andreev reflection in bosonic condensates 
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We study the bosonic analog of Andreev reflection at a normal-superfluid interface where the 
superfluid is a boson condensate. We model the normal region as a zone where nonlinear effects can 
be neglected. Against the background of a decaying condensate, we identify a novel contribution 
to the current of reflected atoms. The group velocity of this Andreev reflected component differs 
from that of the normally reflected one. For a three-dimensional planar or two-dimensional linear 
interface Andreev reflection is neither specular nor conjugate. 
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Andreev reflection occurs at a normal-superconductor 
interface when an incident electron/^hole on the normal 
side is reflected as a hole/electron [Ij. It is a current- 
carrying process which, within a mean-fleld description, 
conserves the quasiparticle but not the charge current. 
The Andreev reflection of an incident hole is equivalent to 
the emission of a singlet electron pair from the supercon- 
ductor into the normal metal [4I . It has been observed in 
superconductors Q and superfluid ^He J| , and has been 
proposed for fermionic quantum gases (5|. In these Fermi 
superfluids the conversion between particles and holes re- 
flects the fermion-pair structure of the superfluid phase. 
Andreev reflection determines the transport properties of 
fermionic superfluid-normal interfaces. It seems natural 
to investigate its bosonic analog. Andreev-like processes 
have been proposed for cold atom realizations of Bose- 
Hubard models [6| and Luttinger liquids [7|. These two 
works investigate the generation of reflected dips in the 
particle density following the incidence of a density bump 
on an interface. 

Here we investigate the possible analog of Andreev re- 
flection at the interface between a bosonic superfluid and 
a normal gas. We seek an analog closer to the fermionic 
case by exploring the conversion between different types 
of Bogoliubov quasiparticles, which for bosons are defined 
through the transformation ip = e"^^yj(u„7n — v:^7n)y 
where ip is the field fluctuation operator, (jj is the conden- 
sate phase, and 7„ destroys quasiparticle n [8| . In bosonic 
flelds the u and v components play asymmetric roles; for 
instance, under usual circumstances ti = in a normal 
gas. In this context, we define Andreev reflection in a 
normal-superfluid (NS) interface as the scattering event 
which generates a quasiparticle with w ^ on the N side 
following the incidence of a conventional {v = 0) quasi- 
particle from the same side. The hole-like quasiparti- 
cle (which however contributes positively to the density) 
travels at a speed different from the normally reflected 
one, which should be experimentally observable. 

A natural way of forming a NS interface in boson gases 
is to create a potential step that forces a difference in the 
atom densities on each side of the step. At a given tem- 
perature, the densities may be such that the gas is normal 



on the left and superfluid on the right. However, since 
in a normal boson gas the chemical potential stays below 
zero kinetic energy, transport near equilibrium is possi- 
ble only if the condensate is confined. The v component 
represents propagation below the chemical potential, so it 
can only be evanescent on the N side within that scenario 
[Sl, llO| . As Andreev reflection from a confined condensate 
is ruled out, one is led to consider the alternative case of 
a decaying condensate. 

We define the normal (or non-superfluid) region as that 
in which the condensate flows faster than the local speed 
of sound; in our case, the outgoing coherent beam. For 
simplicity, we assume the interaction coupling strength 
to be zero on the N side, which does not change the 
essential physics. In general, transport through an in- 
terface separating subsonic from supersonic flow poses a 
new paradigm in superfluid transport |11'|. 

Like for fermions, particle current is not conserved 
within a (non-selfconsistent) mean-fleld description of 
Andreev processes [ij, ll3|. In both cases the conden- 
sate is responsible for the loss or gain of particles. How- 
ever, there are several important differences with respect 
to the fermion case. The eigenvalue problem whose solu- 
tions are the Bogoliubov quasiparticles is non-Hermitian. 
This results in a peculiar normalization condition [v = 
fHitp — |tip)dx = 1] for the quasiparticle wave functions 



[14l | . We will see that the conservation of ly yields uncon- 
ventional relations between the scattering amplitudes. A 
counterintuitive but straightforward consequence is that, 
unlike for fermions, both particle-like and hole- like ex- 
citations carry an increase in the density (with respect 
to the Bogoliubov vacuum), since p ~ JMp -I- |z;p. The 
existence of these differences suggests that bosonic quasi- 
particles define a rich novel class of quantum transport 
problems. 

The need for a confined condensate to generate An- 
dreev reflection can be better appreciated from a study of 
the Bog oliubov - de Gennes (BdG) equations for bosons 
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FIG. 1: Dispersion relation on the normal side for branches 
u (particle-like) and v (hole-like) for a confined (left) and a 
decaying (right) condensate. Labels are explained after Eq. 
(Ul in main text. Elastic conversion between propagating 
branches (Andreev reflection) requires 5fi > 0, i.e. a decaying 
condensate. 
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(1) 



C{x) = -\Dl + V{x) + 2nQ{x)g{x) - /i 



Here, h ~ m ~ 1 units have been used, no(x) is the 
condensate density, g{x) is the (piecewise constant) cou- 
phng strength, Dx = dx — i4>' with (j){x) the condensate 
phase, and /i is the condensate chemical potential. The 
scattering channels are defined by the plane wave solu- 



tions (u, V 
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) in the flat region V{x) 



V^. 



The dispersion relations in the non-superfluid region 
{g{x)nQ{x) — 0) are shown in Fig. [T] In the confined 
case ((5/i = ji — Vo < 0), the left Fig. [1] shows that, 
away from the interface, asymptotic propagation with 
£ > is only possible for the u (particle- like) component, 
while the v (hole-like) component can only be evanes- 
cent. A decaying (5/i > 0) condensate exhibits a richer 
scenario. The right Fig. [1] shows a crossing of the u- 
and w-branches. This permits elastic conversion from a 
particle into a "hole" , in formal analogy with fermionic 
Andreev reflection. For a given energy within the allowed 
crossing range e € {—5fi,Sfi), there are four propagating 
waves. The right (left) going wave numbers are fc^ and 



— fc^ (— fcjv ^^'^ ^n)j where fc^ = y/2{Sii ± e). The 
group velocities are those of a quadratic dispersion rela- 
tion, WjY = fcjy. Thus, for an incident particle from the 
normal side with velocity Wi — •\/2((5/i -t- e), there are two 
outgoing channels with velocity Wf satisfying Wf — Wi 
(normal reflection) or wf + w'i — 4 Sfi (Andreev reflec- 



FIG. 2: Schematic plot of the scattering problem. Against 
the background of the coherent beam emitted by a decay- 
ing condensate, a single atom impinges from the normal side. 
Three scattering processes result: normal reflection, transmis- 
sion into a phonon mode on the superfluid side, and An- 
dreev reflection on the normal side, with group velocities 



respectively. Labels are explained after Eq. 



([T]) in main text. 



tion). The prediction of this novel form of propagation 
on the normal side is the main result of this paper. 

On the superfluid side, excitations travel with wavevec- 
tor kg = \f2\\J e^ -f /i^ — l>^ ■ Here, /i — gnj^, with 
n± = limj;^±oo no(x). At low energies, the speed ws = 
de/dkg approaches the speed of sound c = gn^ . An An- 
dreev transmission process occurs for energies e G (0, 5^) 
when a phonon incident from S is transmitted into N as 
a hole-like {u — 0) quasiparticle. The other transmis- 
sion channel is particle-like [v — 0) and travels in N at 
a different velocity {wj^ vs. w^). Remarkably, the out- 
going hole-like scattering channels which characterize an 
Andreev process may exist despite its anomalous normal- 
ization 1/ < 0. We note, however, that the conservation 
of V applies to the total scattering state: a particular 
channel may have v < Q whereas the total state retains 
the standard normalization u > Q. 



Bogoliubov quasiparticles are generally defined against 
the background of a particular condensate wave function. 
We consider a one-atom potential of the type shown in 
Fig. [21 which provides a model for a leaking conden- 
sate. The condensate wave function ^ satisfies the time- 
dependent Gross-Pitaevskii equation 



V{x) 



Z5{x) + V(ie{-~x) . 



(2) 
(3) 



which is constant for any two stationary solutions of Eq. 
^ with same energy e. When appHed to V'jv and if^s, the 
constancy of ([6|) generates the following relations between 
the scattering amplitudes: 



In the asymptotic regions it admits a solution of the 
type ^(x ^ ±00) ~ g«(?2;-Mt)_ In the previous discussion 
on scattering channels, we have assumed q = Q for sim- 
plicity, which — it can be shown — does not change the 
main conclusions. 

A stationary solution to ^ and ^ can be found 
analytically 17|. An important question is whether or 
not that solution is stable. For Zj ^JJL ^ 1, the exact 
stationary solution resembles closely the wave function 
^(a;) — y/n^9{~x) + y/n+9{x). For this approximate, 
step-like condensate, we have proved that the solutions 
to the BdG equations ([T]) do not have complex eigenval- 
ues. In other words, in the thick barrier limit the exact 
stationary solution resembles a solution which has been 
proven to be stable. In general, the condensed atom den- 
sity on the N side is uniform and equal to n_ = \j\/ ^j2EJi, 
with atom velocity w — —\/2SJl, as determined by sim- 
ple energy considerations. In this thick barrier limit, the 
resulting current is j ~ —fj?^/25ii/gZ^. This solution 
applies to an interface separating a supersonic (normal) 
from a subsonic (superfluid) region [ll|. 

Next we write the general form of the scattering states 
within the approximation of a flowless condensate. For 
an atom incident from the N side with \e\ < 6^ we have 





tpxs+{-::)te^ 


+ (l)^aXAr, X ^ -00 

i^'^s ^, a; ^ 00 

(4) 


where x% = c^'^n^^w^, xs = e*«s"(^)4-/^/iI^ and 


kg = V2[-\/e2 -(- ^2 _|_ ^]i/2 and uq, vq are such that Uq + 


4 = \/e' + m' 


/ £ and uo/vo = 


= ^/l + eV/x2 + e/fi. For 



an incoming phonon from the S side with energy < e < 
Sfi, we have 



■0s 



(o)^» (xn) 



O^aXN^ 



{Z)[xs+rr>xs] + {-::> 



00, 



(5) 
The subindeces stand for normal (a), Andreev (n), 
phonon (p) and evanescent (e), while r and t represent 
reflection and transmission coefficients. Equations (|4]) 
and ([5]) do not include the case e > Sfi, which shows no 
Andreev processes, or that of an incoming channel with 
v < 0, which does not seem realizable unless the beam 
comes from another Bose condensate. 

Conservation laws may be obtained from the general- 
ized Wronskian 



r„p-|r,|2+sgn(e) |ipp = |t„| 



\ta\ 



^ 1 



fn t^ + r tp Ta t^ — V) 



(7) 

The negative signs here can be traced back to the non- 
Hermitian character of the effective Hamiltonian in the 
bosonic BdG equations ([T]). This represents a major dif- 
ference with respect to the fermionic Bogoliubov problem 
and is a source of difficulties in the interpretation of the 
particle- hole transformation for bosons. 

If we assume a step-like profile for the condensate den- 
sity, the scattering amplitudes can be found analytically. 
Here we note some important trends in the low trans- 
parency limit {Z/y/Ji 3> 1): The Andreev reflection am- 
plitude \ra\ scales like Z^^, which indicates that, like for 
superconductors [18[, an Andreev reflection process in- 
volves the transmission of two atoms. This reflects the 
pairing correlations which exist in the quantum deple- 
tion cloud of the condensate and which are expressed 
through the structure of the Bogoliubov transformation 
and the BdG equations [16J. Analogously, the scaling 
\ta\, \tp\ ^ Z~^ suggests that for large Z the correspond- 
ing scattering processes involve the transmission of one 
atom. 

More information on conservation laws can be obtained 
from the continuity equation 



dt\uf + dxju = dt\v\'^ + dxjv* = 2 nQ{x)g{x) lm{uv*) , 

(8) 
where js = lin{s* D^s), which tells us simultaneously 
about the lack of conservation of the atom mass current 



(density 



current j„ -I- jy* ) and the conser- 



vation of the quasiparticle current (density \u\ — \v\ , 
current ju — jv*), the latter being also expressed by the 
constancy of the Wronskian ^ . Both mass density and 
current are defined as their total values in the presence of 
the excitation of wave function (") minus the correspond- 
ing values of the Bogoliubov vacuum, which includes the 
zero-temperature depletion cloud [8[ and thus a fluctu- 
ating atom flux. We note that, while the u component 
carries the current in the usual way, the v does it in the 
opposite one. Another important feature is that, unlike 
for its fermionic counterpart [l8|, the mass current has 
the same sign as the group velocity in all channels (the 
mass density is always positive). 



W,, 



UiDlu* - u*DxUt + ViDxV* - v*D*u, 



(6) 



As in the case of superconductivity [13|, ll8|, we in- 
terpret the non-conservation of atom number as the ex- 
change of atoms with the Bogoliubov vacuum. The net 
balance of extracted atoms per unit time can be calcu- 
lated by substracting, from the current at a point deep 
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FIG. 3: Andreev reflection at a bosonic three-dimensional 
NS interface. Indeces i, n, a stand for incoming, Andreev re- 
flected and normaUy reflected channels. The arrows represent 
group velocities. The angle of the Andreev reflection channel 
is sensitive to the sign of the energy e. 



in the superfluid region, the current at any point in the 
normal region, A/ = j{x+) — j{x-). We obtain 



A/ = 



'■ + \raf 



\tp? 



2 



-1 



-2 sgn(£) k+ k+ 1 



1 + ey^i^ z2 ■ 

(9) 
For the second equality, which gives the current to lead- 
ing order in Z~^, we have used JTl) and the exact results 
for the scattering amplitudes [17|. Thus, for an atom 
incident from the N side, the scattering process which 
follows is such that atoms are lost to the Bogoliubov vac- 
uum if £ > 0, while they are extracted from that vacuum 
if £ < 0. From the continuity equation ([8]) it is clear that 
the transfer of atoms between the condensate and the 
quasiparticle field takes place exclusively on the S side, 
since g{x) — in N. 

The generalization of the above discussion to the case 
of a three-dimensional planar or two-dimensional lin- 
ear interface yields qualitative differences with respect 
to the superconducting case. Invoking the conserva- 
tion of energy and parallel momentum, we find that for 
£ e {—S^, SfJ.) the beam angles (see Fig. [3]) satisfy the 
relation 

(1 - e/Sfi) tan^ a.j - (1 + e/S^i) tan^ «„ = 2 e/S^i . (10) 

This result leads to two conclusions: (a) If < £ < 6fi, 
no Andreev reflection occurs for tanai < y^e/2{5fi — e) 
and one has aa < ai. (b) If — (5/i < £ < Andreev 
reflection is possible, but only satisfying Ua > on and 
tantta > v^— £/2((S/i-f £). 



The observation of bosonic Andreev reflection poses a 
challenge. A tunable coupling constant can be achieved 
through the handling of Feshbach resonances [19| or, in 
the ID limit, by increasing the width of the conflning 
channel. Such long channels can form in the vicinity 
of a planar semiconductor chip. A potential barrier at 
the interface can be introduced with a properly focused, 
blue-detuned laser. The crucial measurement of escaping 
velocities may rely on precise Bragg spectroscopy or on 
the interference between incoming and outgoing beams. 
The 2D/3D case presents greater potential flexibility be- 
cause anomalous reflection can occur in a range of differ- 
ent directions. Conflncmcnt by a linear interface in 2D 
could be achieved with a laser beam. For instance, ®^Rb 
atoms may be confined in a waveguide with transverse 
trapping frequency of 500 Hz and experience a delta- 
barrier oiZ — hyfii, if they are exposed to a narrow Gaus- 
sian laser of 1 /im waist, blue-detuned by 10 nm to the 
52Si/2-52Pi/2 transition 20| (780.24 nm, linewidth 127r 
Hz, saturation intensity 1.64 mW cm"^). With a step 
barrier such that (5/i = 0.5/i, bulk superfiuid linear den- 
sities in the range 0.2-20 ^vcT^ place the system in the 
low-density mean-field region, where the present estimate 
applies. For that density range we obtain ID coherence- 
length 7-0.7 /xm, laser intensity 80-800 mW cm"^, and 
10-1000 ^m s~-^ for the Andreev velocity scale (^ 2-v/5/i)i 
close to the superfluid sound speed of 70-700 /xm s^^ . For 
these parameters the emitted condensate linear density 
is 0.04 times the bulk S density, while the corresponding 
condensate escape velocity is the same as the speed of 
sound in S, the equality being due to our choice of (5/i. 
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